Impact of dimerization and stretching on the transport properties of molybdenum 

atomic wires 
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We study the electrical and transport properties of monoatomic Mo wires with different structural 
characteristics. We consider first periodic wires with inter-atomic distances ranging between the 
dimerized wire to that formed by equidistant atoms. We find that the dimerized case has a gap in 
the electronic structure which makes it insulating, as opposed to the equidistant or near-equidistant 
cases which are metallic. We also simulate two conducting one-dimensional Mo electrodes separated 
by a scattering region which contains a number of dimers between 1 and 6. The I — V characteristics 
strongly depend on the number of dimers and vary from ohmic to tunneling, with the presence of 
different gaps. We also find that stretched chains are ferromagnetic. 
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I. INTRODUCTION 

The field of Nanoelectronics has emerged as one of 
the most important branches of Nanoscience and Nan- 
otechnology, due to both the fundamental quantum phe- 
nomena occurring at this scale and to its technological 
relevance in the context of miniaturization and perfor- 
mance improvement of electronic devices. Remarkable 
progress in the experimental techniques, achieved during 
the last years, for growing, manipulating and measuring 
at the nanoscale has allowed to demonstrate several pro- 
totype devices such as conducting wires, point contacts 
and switches.— 

Metallic atomic-scale nanowires, in particular, have at- 
tracted considerable efforts devoted to their production 
and isolation^— Also infinite monoatomic wires have 
been theoretically investigated recently in order to bet- 
ter understand the magnetic and transport properties of 
one-dimensional atomic contacts^— ^From such ideal- 
ized models we can gain valuable information which can 
be compared later with the realistic systems to get in- 
sight in the influence of defects, impurities, host material 
or the structure of the contacts. 

Understanding the transport properties like the cur- 
rent flow through these nanoscale devices, has been 
a challenge from the theoretical side due to the non- 
equilibrium quantum kinetic description required. The 
modern theory of quantum transport is based on the 
scattering theory or on the Keldysh-Kadanoff-Baym non- 
equilibrium Green's function (NEGF) formalism. 8 The 
SMEAGOL code2J£ is a flexible and efficient implemen- 
tation of the NEGF formalism. SMEAGOL obtains the 
Hamiltonian from the density functional theory (DFT) 11 
code SIESTAfi 2 - which uses pseudopotentials and a lo- 
calized basis set of pseudo- atomic orbitals, and calculates 
self-consistently the density matrix, the transmission and 
the current for each bias voltage. 

Very recently, individual molybdenum chains have 
been produced and controlled by encapsulating them 
inside carbon nanotubes^ Mo is indeed an interest- 



ing metal element as regards its ability to form one- 
dimensional (ID) structures, even in free environment. 
Recent ab-initio calculations of free-standing Mo atomic 
clusters have shown that linear atomic chains containing 
up to four atoms are considerably more stable than two- 
and three- dimensional structures^. This result is rele- 
vant in the context of Mechanically Controllable Break 
Junction experiments (MCBJE), where atomic chains 
are fabricated at ultra-low temperatures by bending and 
eventually breaking up an atomic strip. Interestingly, 
only the elements gold, platinum and iridium have been 
found to make thenijiS, although other elements like Sn 
could be possible candidates?^. Due to its tendency to 
make linear clusters^ Molybdenum could also be used 
to make atomic chains in MCBJE. Therefore, a trans- 
port calculation of such chains is indeed timely. 

Besides, as a consequence of Mo having an exact half 
band filling, the geometry is characterized by the forma- 
tion of tightly bound dimers, which lead to a insulating 
ground state. Further, these chains show a Peierls-type 
insulator-metal which is activated when the two inter- 
atomic distances are varied so that the dimerization is 
lost. This article analyzes the impact of this metal insu- 
lator transition in the transport properties of the chains. 

More specifically, we study the electrical and transport 
properties of infinite monoatomic Mo wires with differ- 
ent structural characteristics. We consider first periodic 
wires with inter-atomic distances ranging between the 
dimerized wire to that formed by equidistant atoms. We 
have also simulated two conducting one-dimensional Mo 
electrodes separated by a scattering region which con- 
tains a number of dimers between 1 and 6. Since the 
conductance is closely related to the electronic struc- 
ture, which depends dramatically on the distance be- 
tween atoms, interesting effects are expected to appear 
as a function of the structural changes induced by the 
dimerization. 

In the next section we give the details of our DFT cal- 
culations within the SMEAGOL code. In Section III we 
present the results obtained for the infinite monoatomic 



Mo wires between the dimerized and the equidistant con- 
figurations. Section IV is devoted to the scattering region 
formed by finite dimerized chains anchored to infinite 
wires with equidistant atoms. In section V we discuss 
the results for a Mo wire with interatomic distances sim- 
ilar to those experimentally observed when the wire is 
encapsulated in the nanotubes. Magnetic effects are ex- 
pected here due to the considerable expansion of the in- 
teratomic distances. Finally, Section VI summarizes the 
main conclusions. 



II. DETAILS OF THE DFT APPROACH 

We calculated the exchange and correlation potential 
with the generalized gradient approximation (GGA) as 
parametrized by Perdew, Burke and Ernzerhof,— We re- 
placed the atomic Mo core by a nonlocal norm-conserving 
Troullier-Martins 1 - pseudopotential which was factorized 
in the Kleinman-Bilander formic and generated using the 
atomic configuration 4d 5 5s 1 5p°, with cutoff radii of 1.67, 
2.30 and 2.46 a.u., respectively. We also included nonlin- 
ear core corrections, generated with a radius of 1.2 a.u, to 
account for the significant overlap of the core charge with 
the valence d orbitals and avoid spikes which often ap- 
pear close to the nucleus when the GGA approximation is 
used. We tested that this pseudopotential reproduced ac- 
curately the eigenvalues of different excited states of the 
isolated Mo atom. We employed a linear combination of 
pseudoatomic orbitals to describe the valence states of 
our chains ^ The basis set included double- C polarized 
(DZP) orbitals, i.e. two radial functions to describe the 
5s shell and another two for the Ad shell, plus a single 
radial function for the empty p shell. We used a numer- 
ical grid defined by an energy cutoff of 250 Rydberg to 
compute the exchange and correlation potential, and to 
perform the real-space integrals that yield the Hamilto- 
nian and overlap matrix elements. We also smoothed the 
Fermi distribution function that enters in the calculation 
of the density matrix with an electronic temperature of 
300 K and used a conjugate gradients algorithm, 21 to re- 
lax the atomic positions until the interatomic forces were 
smaller than 0.01 eV/A. Finally, we performed careful 
tests for particular cases to ensure the quality and sta- 
bility of the basis set and the real space energy cutoff 
employed. We found that the results were hardly modi- 
fied when the DZP basis was replaced by a triple-^ dou- 
bly polarized basis. Similar results were also obtained by 
considering an electronic temperature of 100 K. 



III. INFINITE PERIODIC MONOATOMIC Mo 
WIRES 
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The ground state of the periodic monoatomic Mo wire 
is formed by tightly bonded dimers of 1.58 A inter-atomic 
distance and 3.01 A inter-dimer distance. The resulting 
wire is non-magnetic, since the Mo atom has exactly a 



FIG. 1: Transmission as a function of the energy referred to 
the Fermi level for a Mo atomic wire with interatomic dis- 
tances of (a) 1.58 A and 3.01 A (b) 1.74 A and 2.85 A (c) 
1.88 A and 2.71 A (d) 2.02 A and 2.57 A (e) 2.16 A and 2.43 
A (f) equidistant 2.30 A 
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FIG. 2: Density of States projected on the hybrid cartesian 
s-, p- (first column) and d- (second column) orbitals, as a 
function of the energy referred to the Fermi level for the same 
wire configurations as those of Fig. ([TJ). 



half filled d-shcll, which is a consequence of the strong 
covalent bond formed between the two Mo atoms of the 
dimer to achieve a closed-shell electronic configuration. 
Notice that dimer formation arises in finite length chains 
containing as few as 4, 6 and 8 atoms, where one can not 
strictly speak of electronic bands. This phenomenon is 
therefore a Peierls dimerization transition driven by the 
electron-phonon coupling.— Indeed, the magnetic nature 
of these chains is paramagnetic, as opposed to the antifcr- 
romagnetic character which would result by using nesting 
arguments. Further, chains with an odd number of Mo 
atoms (3, 5, 7) also dimerize, the dimers being in a singlet 
state, and the unpaired atom providing 5/is. 

These chains are electronic insulators. To further qual- 
ify their insulating behavior, we have used SMEAGOL to 
compute the transmission coefficients. SMEAGOL com- 
putes the current using the Landauer formula 8 
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FIG. 3: (color online) Transmission (first column) and real 
space plot of the charge density around the Fermi level (sec- 
ond column) for a scattering region with a number of dimers 
ranging from 1 to 6. 



I(V) = Y J *ET(E,V) (f L (E,V)- f R (E,V)) (1) 

where T(E, V) are the energy- and voltage-dependent 
transmission coefficients of the junction, and 
Jl/r{E,V) = f(E — /xl/r) are the Fermi functions 
of the left/right electrodes, whose chemical potentials 
A*l/r = A 4 i eV/2 are the equilibrium chemical potential 
shifted by the voltage bias V. The low-energy and 
low-voltage transmission coefficients can be expanded as 



T(E,V) ~T ,+T E E + T V V + 0(E 2 ,V 2 ) (2) 

where T = T(0,0), and T Ey = g^yj(E = 0, V = 0). 
As a consequence, the low-voltage differential conduc- 
tance is equal to 



G{V) 
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G (T(0, 0) + 27V V) ca G T(0, 0) (3) 




Number of dimers in the scattering region 

FIG. 4: Logarithm of the transmission at the Fermi level as 
a function of the number of dimers in the scattering region. 

where Go = 2e 2 /h is the conductance quantum unit. 
Notice that the first-order energy term in the Taylor ex- 
pansion drops because of the symmetry properties of the 
integral. Therefore, the zero-energy, zero- voltage trans- 
mission coefficients provide an estimate of the low- voltage 
differential conductance. Furthermore, for symmetric 
junctions such as those that will be studied in this ar- 
ticle T(E, V) = T(E, -V), so that G(V) = G(-V). The 
differential conductance must be flat at very low voltages, 
and only second-order voltage terms contribute. 

We plot in Fig. ([T]) (a) the zero-bias transmission 
as a function of energy T{E, 0) for the dimerized infi- 
nite Mo atomic wire. This chain does not have conduc- 
tion channels in a window of about 1.5 eV around the 
Fermi energy, which is confirmed by plotting the Density 
of States (DOS) N[E). Notice that in a ballistic one- 
dimensional conduction model the zero- voltage transmis- 
sion coeficient at a given energy can be written as 

T(E)(xJ2 v n(E)N n (E) (4) 

n 

where the sum is over all the electron bands, and v n {E) 
and N n (E) are the group velocity and the density of 
states at energy E of each band. Since in one dimen- 
sion, N n = l/(Tiv n ), the transmission coefficient T(E) 
just counts the number of bands at energy E, which are 
called transmission channels. For energies below -0.5 eV 
the transmission exhibits a plateau with a magnitude of 
1, which corresponds to a single low-scattering sp band, 
or channel, and for energies above 1 eV, T = 2 due to 
contributions from the antibonding d orbitals. 

In order to understand the interplay between dimeriza- 
tion and transport, we have calculated the zero- voltage 
transmission coefficients T(E) for six chains with differ- 
ent inter-atomic distances. The unit cell of all chains 
contains two atoms, and has a lattice constant a = 4.6 
A. We have varied the inter-atomic distance (d) of the 
two atoms in the unit cell from d\ = 1.58 A (ground 
state dimcr) to d-2 = 2.30 A (equidistant atoms). To 
characterize all these chains, we introduce the dimer or- 
der parameter A = 2d/a, which varies from A = 0.69 for 
d = d-y to A = 1 for d = e?2 • 
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FIG. 5: I — V characteristics for a chain with a scattering 
region with a number of dimers ranging from 1 to 6. 
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FIG. 6: Transmission (a) and projected density of states on 
the hybrid cartesian s-, p- (b) and d- (c) orbitals, for a Mo 
chain with an interatomic distance of 3.5 A. The energies are 
referred to the Fermi level. 



We note first that the ground state of all chains is non- 
magnetic. Figs, (p} (a) to (f) illustrate the evolution of 
the transmission coefficients with the order parameter. 
What we find is that as A increases from 0.69, many oc- 
cupied and unoccupied channels move gradually towards 
the Fermi level, increasing the conductive nature of the 
chains. There is an insulator-metal transition at about 
A = 0.88 (Fig ([I]) (d)) and for larger order parameters 
the chains become metallic. In other words, the insu- 
lating character is retained for dimerized chains, where 
the dimer is strained [A S (0.69,0.88)]. Conversely, not 
only perfect (A = 1), but also slightly dimerized chains 
[A G (0.88, 1)] are metallic. Notice also that the width of 
the electronic bands and the plateaus of T(E) increases 
with A, indicating that the electronic states become more 
delocalized. The actual conductance G — T(0) Gq in- 
creases step- wise from Go at A = 0.88 to 6Gq at A = 1. 

To understand how each state participates in the con- 
ductance, we plot the evolution of the projected density 
of states for each s-, p- and d-orbital in Fig. ([2]). Notice 
the chains are oriented along the z-axis, which is also the 
transport direction in our simulations. The d 32 2_ r 2 or- 
bitals have therefore the largest hybridizations (that is to 
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FIG. 7: Transmission and projected density of states on the 
hybrid cartesian s-, p- and d- orbitals for the up (a) (c) (e) 
and down (b) (d) (f) spin contributions, for a magnetic Mo 
chain with an interatomic distance of 3.5 A. The energies are 
referred to the Fermi level. 



6 



say, Hamiltonian matrix elements) among all d-orbitals, 
followed by the d zx and d zy orbitals (which are equivalent 
by symmetry), leaving the d x a_ y a and d xy as the less hy- 
bridized orbitals. As a consequence, the d 3z 2_ r 2 orbital 
feels more strongly the dimerization of the chain, while 
those states related to the d x 2_ y 2 and d xy orbitals are 
the least affected. The most relevant effect of the dimer- 
ization is the opening of a gap between the bonding and 
antibonding states, which is largest for the d3 Z 2_ r 2 or- 
bital (6 eV for A = 0.69) and smallest for the d x 2_ y 2 
and d xy orbitals (3.5 eV for A = 0.69), while the gap is 
of 5 eV for the states originated from the d xz and d yz 
orbitals. There is also a hybridized sp z orbital below 
the Fermi level, and a p x and p y contribution hybridized 
with the d xz and d yz orbitals. When the order parame- 
ter A approaches 1, the first channel which reaches the 
Fermi level and participates in the zero-bias conductance 
is the sp z , followed immediately by the d xz and d yz . In 
contrast, the gap belonging to the d 3z 2_ r 2 orbital only 
closes when A ~ 1, and therefore only contributes to the 
zero-bias conductance when the atoms in the chain are al- 
most equidistant. When the chain is completely equidis- 
tant the gaps in the densities of states disappear and 
all the d-orbitals contribute to the conductance. Thus, 
the 6 channels of conductance at the Fermi level for the 
equidistant chain comes from the 5 d-orbitals and from 
the sp z -orbital, since the contribution of the p^ and p y 
in the d-orbitals is negligible. 



IV. ONE-DIMENSIONAL Mo ELECTRODES 
SEPARATED BY A DIMERIZED REGION 

In order to investigate the role of dimers in the scat- 
tering region between two conductive electrodes we have 
considered the formation of 1 to 6 dimers within an infi- 
nite equidistant wire. Thus, at both sides of the scatter- 
ing region we have ideal conductive ID electrodes. The 
question that arises is if the current flow is disrupted 
and how by the dimers. The scattering region, which we 
call the extended molecule, includes also part of the elec- 
trodes (four atoms at each side) to ensure the continuity 
of the potential with that of the infinite electrodes. To 
determine the I—V curve we need to evaluate the charge 
density by using the lesser Green's function (instead of 
the retarded one, since electrodes at different potential 
establish a nonequilibrium population in the extended 
molecule) and the self-energies of the electrodes at differ- 
ent potentials for all the energy range in the bias window 
(not only at E-p ) . The intensity of the current at a given 
bias is obtained by using Eq. (1), where 

T(E, V) = tv[T L G^T R G^](E, V). (5) 

= [e+Su - H M - - Eg]" 1 is the retarded 
Green's function of the extended molecule, with e + = 
lim,5^o+ E + iS and E^dO ^ ne retarded self-energy of the 



left (right) electrode, and ^l(r) — *Pr(l) — ^l(r)J * s 
the Gamma matrix, which gives the strength of the cou- 
pling between the scattering region and the left (right) 
electrode. 

We plot in Fig. (J3J the zero-voltage transmission for 
different scattering regions with an increasing number of 
dimers, from I to 6. Each panel has attached to the right 
a real-space plot of the density of states integrated in a 
window of 0.5 eV about the Fermi energy. We find that, 
even if we introduce only one dimer, the transmission 
at the Fermi level falls from 6 to only I Go, due to the 
lack of transmission channels in the scattering region. 
For only one dimer we find that some electronic density 
spreads around the dimer, but becomes rapidly localized 
when we increase the number of dimers. This indicates 
again that the transport must be in the tunneling regime, 
where the low bias conductance decreases exponentially 
as a function of the distance between electrodes, since 
electrons at the Fermi energy can only hop from dimer 
to dimer via tunneling events. The transmission T{E) 
begins to resemble that obtained in the case of the infinite 
dimerized wire (Fig. (JTJ , top panel) when the number of 
dimers increases. 

To confirm the tunneling behavior of the conductance, 
we plot in Fig. @ the logarithm of the zero-voltage 
transmission at the Fermi level in a logarithmic scale as a 
function of the number of dimers in the scattering region. 
The points are indeed fitted by a linear function, so the 
transmission coefficients fit to an exponential. To find 
out the effective height V e g of the tunneling barrier, we 
set 



T(0,0) .ieV^^ = 1.224e- 110JV (6) 

where d = 1.58 A is the dimer length, m is the electron 
mass and TV is the number of dimers in the chain. Notice 
that neither m nor V e ff are properly defined separately. 
Rather, only their product m x V e tf can be extracted 
from the fitting. We pick a common practice when fitting 
conductance curves in vaccum tunneling, and select m to 
be the physical mass of the electron, thus dumping all 
the non-trivial renormalizations into the effective height 
of the barrier. We find that V e g ~ 1.8 eV, which agrees 
reasonably well with the gap in the DOS of an infinite 
dimerized chain. 

In Fig. §5§ we plot the I—V characteristics in a voltage 
range between -2 and +2 volts, for the different scatter- 
ing regions mentioned above. We note that the I — V 
curves look like featureless Ohmic curves when the scat- 
tering region comprises just one or two dimers. For a 
larger number of dimers the I — V curves start to de- 
velop a gap of about 2 volts, and a number of additional 
high-voltage shoulders. These evolve into negative dif- 
ferential resistance (NDR) characteristics for large num- 
bers of dimers. The NDR is produced by a small peak 
in the transmission coefficients T(E 1 V), which is located 
around -0.75 eV. This peak appears for high absolute val- 
ues of V, but its height has a sudden drop when the V is 
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increased beyond 1.7 volts. The drop leads to a concomi- 
tant decreases in the electrical current. We find that if 
V is increased even further, a second peak enters inside 
the integration window when V ~ 1.9 volts, so that the 
electrical current increases again. The emergence and 
evanescence of transmission peaks can be understood in 
terms of the disorder introduced by the bias-dependent 
on-site energies^. We have found a similar NDR effect 
in thiol-capped polyynes junctions. 23 

V. STRETCHED MONOATOMIC Mo WIRES 

The interatomic distances d recently measured be- 
tween Mo atoms in individual monoatomic chains en- 
capsulated inside carbon nanotubesi 3 - ranged from 3.2 
A to 3.8 A. Those wires are then far from being dimer- 
ized; rather, they are stretched, since those distances are 
much longer than 2.3 A, which corresponds to equidistant 
free-standing chains. 

We expect that isolated chains with such long d become 
ferromagnetic because of the Stoner criterion: they must 
have narrower bands and, as a consequence, a high den- 
sity of states at the Fermi level in the non-magnetic state. 
Indeed, we have found that the ground state configura- 
tion is characterized by a full spin polarization, where the 
magnetic moment is saturated to its maximum possible 
value of 6 fib per atom, and by a large exchange splitting 
of xx eV. Hence, all 5d as well as the 6s orbital are com- 
pletely spin polarized. This ferromagnetic ground state 
is more stable than the paramagnetic state by 2.3 eV. 

We expect that this spin polarization must have an im- 
portant impact in the transport behavior. Indeed, Hafner 
et al2^ found in a recent work that atomic contacts of 
Fe, Co and Ni carry spin polarized currents, where the 
degree of the spin-polarization of the current depended 
on the type of transport regime (contact or tunneling. To 
estimate the extent of this impact, we have calculated the 
transmission in both the non-magnetic and the magnetic 
states. 

We have computed first the transmission coefficients 
of an isolated Mo wire in the non-magnetic state with 
equidistant interatomic distances and a unit cell similar 
to that found experimentally. The results are shown in 
Fig. ©. We find that all six s- and d-channels con- 
tribute to the transmission coefficient at the Fermi level, 
like in the equidistant regime discussed in section III. 
This might look surprising at first glance since the dis- 
tance between atoms in this stretched chain is even larger 
than the inter-dimer distance of the dimerized chains 
with insulating behavior discussed in section III. How- 
ever, we note that in the equidistant stretched chain 
the d-orbitals do not form covalent bonds, and since no 



bonding and antibonding states are formed (with their 
corresponding gap), the d-electrons are free to partici- 
pate in the transmission at the Fermi level. Away from 
the Fermi level, however, the transmission coefficient de- 
creases much faster than in the previous case, since weak 
bonds due to electron localization generate small band- 
widths and, therefore, small energy intervals in the trans- 
mission. The narrowing of the bandwidth due to electron 
localization is clear if we compare the densities of states 
of Fig. ((2} with those corresponding to the equidistant 
wire in Fig. (JBJ (lower panel). 

We find that the ferromagnetic state is insulating, cur- 
rent is fully spin-polarized for the ferromagnetic state, 
since all s and d channels with spin down (minority elec- 
trons) between -3 and 3 eV contribute to the transmis- 
sion above to the Fermi level, as can be seen in Fig. (J7J. 
We do not observe half-metallic behavior at the Fermi 
level, however, due to the presence of the large spin split- 
ting which opens to completely occupy the d and s levels 
with unpaired electrons. The encapsulation inside a nan- 
otube is expected to modify these electronic and trans- 
port properties, since there will be a net charge transfer 
between the chain and the nanotube. 



VI. CONCLUSIONS 

We have found that dimerization of molybdenum 
atomic wires has dramatic effects on their electronic and 
transport properties. While equidistant wires are metal- 
lic and have a very high zero-bias conductance, dimerized 
wires show a large gap which makes them insulating. As a 
consequence, the transport properties of dimers between 
equidistant conducting wires show tunneling behavior. 
The conductance decreases exponentially as the number 
of dimers increases in the scattering region and the I — V 
characteristics start developing a large gap around zero 
voltages. Other features like plateaus and negative differ- 
ential resistance also appear for large numbers of dimers. 

Equidistant stretched chains are highly conductive in 
the non-magnetic case. However, the most stable con- 
figuration is ferromagnetic, with a very high magnetic 
moment of 6 fiB per atom. 
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